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CN ■ Abstract 
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We show that, within a meson-exchange dynamical model describing most of 
the existing pion electromagnetic production data up to the second resonance 

o 



region, one is also able to obtain a good agreement with the 7r° photo- and 

electroproduction data near threshold. The potentials used in the model are 

O . 

derived from an effective chiral Lagrangian. The only sizable discrepancy be- 



X 



tween our results and the data is in the P— wave amplitude P3 = 2M\ + + Mi_ 
where our prediction underestimate the data by about 20%. In the case of 
7T° production, the effects of final state interaction in the threshold region are 
nearly saturated by single charge exchange rescattering. This indicates that 
in ChPT it might be sufficient to carry out the calculation just up to one-loop 
diagrams for threshold neutral pion production. 

Chiral perturbation theory (ChPT) provides us with a systematic scheme to describe 
the low energy interactions of Goldstone bosons among themselves and with other hadrons, 
because it is based on a low energy effective field theory respecting the symmetries of QCD, 
in particular chiral symmetry. There is generally good agreement between the ChPT pre- 



dictions and experiments [T|]. ne case which has been very intensively studied is 7r° electro- 
magnetic production of neutral pions near threshold where very precise measurements @-|8j 
have been performed and the ChPT calculation to one loop 0(p 3 ) (0(p A ) in the case of pho- 
toproduction) has been carried out in the heavy baryon formulation PJTOfl. Nice agreement 
between theory and experiment was reached not only for the S— wave multipoles E 0+ and 



Lq + but for the P— wave amplitudes [pl,|10|l as well. 

As in ChPT, meson-exchange models also start from an effective chiral Lagrangian. 
However, they differ from ChPT in the approach to calculate the scattering amplitudes. In 
ChPT, crossing symmetry is maintained in the perturbative field-theoretic calculation, and 
the agreement with low energy theorems and the data is to be expected as long as the series 
is well convergent. In meson-exchange models, the effective Lagrangian is used to construct 
a potential for use in the scattering equation. The solutions of the scattering equation will 
include rescattering effects to all orders and thereby unitarity is ensured, while crossing 



symmetry is violated. Such models [|Il|-[17j have been able to provide a good description of 
ttN scattering lengths and phase shifts in S, P, and D waves up to 600 MeV pion laboratory 
kinetic energy. 

Meson-exchange models have been constructed for pion electromagnetic production as 



well J15| , |llt i21| and good agreement with the data has also been achieved up to 1300 MeV 
total ttN cm. energy. However, the predictive power of the meson-exchange model for 
electromagnetic pion production near threshold has not been fully explored even though 
the importance of final state interaction (FSI) for threshold tt° photoproduction had been 



demonstrated in several dynamical model studies |12|,[^-|24| prior to one-loop calculations 

of ChPT m. 



In this paper we will use the dynamical model (DM) recently developed in Ref. |25| where 
the dominant FSI effects in the nonresonant contributions at threshold are evaluated using 
the 7rN meson-exchange model developed in Ref. ||13| . Contributions which are related to 



the excitation of resonances are considered phenomenologically using standard Breit-Wigner 



forms. In our previous work |25| such an approach gives an excellent description of pion 



photo- and electroproduction in the first resonance region. Here we apply this model to these 
reactions in the threshold region and compare our predictions with the recent experimental 
data 0-§| for the S— and P— wave multipoles and cross sections, and with the result of 
ChPT. 

In the dynamical model for pion photo- and electroproduction [^,^], the t-matrix is 



given as 



t yn (E) = v yn + v yn g (E) t nN (E) 



(1) 



where v 77r is the 77? transition potential, go and t n N are the nN free propagator and t— matrix, 
respectively, and E is the total energy in the cm. frame. In the present study, the matrix 
elements t n N are obtained in a meson-exchange ttN model [13[ constructed in the Bethe- 
Salpeter formalism and solved within Cooper- Jennings reduction scheme [23]. Both v n N 



and t> 77r are derived from an effective Lagrangian containing Born terms as well as p— and 



lu— exchange in the t— channel p9|j30|j . For pion electroproduction we restore gauge invariance 
by the substitution, 

k-J 



T — > T — h 

<Jfj, ^ <Jn V L2 ' 



(2) 



where J M is the electromagnetic current corresponding to the background contribution of 



J^jTX • 



For the physical multipoles in channel a = {l,j}, Eq. (U) gives [[27 

,Ra(qE,q')v a (q',k) 



t a (<lE, k) = exp (i8 a ) cos5 a 



v a (qE,k) + P dq 
Jo 



(3) 



E(q E ) - E(q') 

where 8 a and R a are the irN phase shift and reaction matrix, in channel a, respectively, q E 
is the pion on-shell momentum and k —\ k | the photon momentum. In order to ensure the 
convergence of the principal value integral, we introduce a dipole-like off-shell form factor 
characterizing the finite range aspect of the potential, f(q,qE) = (A 2 + g| 7 ) 2 /(A 2 + g 2 ) . 
The value for the cut-off parameter, A = 440 MeV, was obtained in our previous work from 
an analysis of the A(1232) resonance multipole M 1+ over a wide energy range. 



For 7T photoproduction, we first calculate the multipole Eq + near threshold by solving 
the coupled channels equation within a basis with physical pion and nucleon masses. The 
coupled channels equation leads to the following expression for the pion photoproduction 
t-matrix in the 7T°p channel: 

t 17T o(E) = Vj„°(E) + v.y n o(E) g 7T o p (E) t 7r o p _ +7r o p (^) 

+ v 77r + (E) g n + n (E) t n+n ^ n0p (E) , (4) 

where t^Op^^Op and t n + n ^ n o p are the nN t-matrices in the elastic and charge exchange chan- 
nels, respectively. They are obtained by solving the coupled channels equation for ttN 
scattering using the meson-exchange model [fl3 ] . Results for ReE 0+ obtained from Eq. (|j) 



are given in Fig. 1, where the FSI contributions from the elastic and charge exchange chan- 
nels (second and third term in Eq. (|4])), are shown by the short-dashed and dash-dotted 
curves, respectively, while the dotted curve corresponds to the LET results, without the 
inclusion of FSI. Our results clearly indicate that practically all of the final state interaction 
effects originate from the ir + n channel. Note that the main contribution stems from the 
principal value integral of Eq. (f|) . 

In the coupled channels approach considered above, the t n N matrix contains the effect 
of ttN rescattering to all orders. However, we have indeed found that only the first order 
rescattering contribution, i.e. the one-loop diagram, is important. This result is obtained by 
replacing the ttN scattering t-matrix in Eq. (|]) with the corresponding potential v^n. As 
can be seen in Fig. 1, the thus obtained results given by the long-dashed curve, differ from 
the full calculation (solid curve) by 5% only. This indicates that the one-loop calculation in 
ChPT could be a reliable approximation for 7r° production in the threshold region. 

If the FSI effects are evaluated with the assumption of isospin symmetry (IS), i.e., with 
averaged masses in the free pion-nucleon propagator, the energy dependence in ReE 0+ in 
the threshold region is very smooth. Below tt + threshold the strong energy dependence (cusp 
effect) [f7||?T| only appears because of the pion mass difference and, as we have seen above, 



is related to the coupling with the ir + n channel. In most calculations, the effects from the 



pion mass difference below the 7r + production channel are taken into account by using the 
K- matrix approach [|I~0|,|32]1 , 



Re Eg = ReEg(IS)-a 7rN uj c ReEg (IS) 



\ 



^ (5) 



where uo and u c are the ir° and n + cm. energies corresponding to E = E p + E 7 and m n + 
m. K +, respectively, and a n N = 0.124/771^+ is the pion charge exchange threshold amplitude. 
Eq1_ (IS) is the 7r°' + photoproduction amplitude obtained with the assumption of isospin 
symmetry (IS), i.e., without the pion mass difference in Eq. @. Such an approximation 
is often used in the data analysis in order to parametrize the Pq+ multipole below 7i + n 



threshold in the form of E 0+ (E) = a + bJl — (uj/uj c ) 2 . In Fig. 2 the results obtained 
from such an approximation scheme are represented by the solid curve and compared to the 
exact coupled channels (dashed curve) and ChPT calculations (dash-dotted curve) 0. We 
see that above n + threshold the difference between the results obtained with physical and 



isospin bases is very small. This is consistent with the finding of Ref . |32| . At lower energies 
the difference becomes visible only very close to 7r°p threshold where the two approaches 
differ by about 10%. In general we can conclude that the simple expression of Eq. (|5|) is 
a good approximation for the pion- mass difference effect, and in the following calculations 
we will use this option to analyze the experimental data. Note that the correct threshold 
dependence of the imaginary part can be obtained from the Fermi- Watson theorem if in the 
threshold region the irN phase shift is taken as a linear function of the ir + momentum, i.e., 
vanishing as q n + — > 0. 

In Fig. 3, we compare the predictions of our model for the differential cross section with 
recent photoproduction data from Mainz 0||. The dotted and solid curves are obtained 
without and with FSI effects, respectively. It is seen that both off-shell pion rescattering 
and cusp effects substantially improve the agreement with the data. This indicates that 
our model gives reliable predictions also for the threshold behaviour of the P-waves without 
any additional arbitrary parameters. As an example, numerical values for E + (in units 
of 10 _3 /m^+) and P-wave multipoles (in units of 10 _3 g/m 3 +) at ir° threshold are given in 



Table I. For the P-wave multipoles we give values for the following linear combinations: 
P 1 = 3E 1+ + M 1+ - Mi_, P 2 = 3E 1+ - M 1+ + Mi_ and P 3 = 2M 1+ + M x _. 

The contributions of Born terms, vector meson exchange, FSI and resonances 
(Sii(1535), A(1232) and Pn(1440)) are listed in Table I to show their relative importance. 
One observes that p— and mostly u— exchange give important contributions to the P-wave 
amplitudes, especially to P3, where Born terms contribute only 2%. Half of P3 comes from 
the vector meson exchanges and the rest from FSI and resonance contributions. Due to the 
large contribution from these three mechanisms, P3 becomes comparable to Pi and Pi- Very 
close to threshold, our model predicts | P3/P2 |— 0.9, and a small but negative value for 
the photon asymmetry £(#„■ = 90°) ~| P3 | 2 — | P2 | 2 at fixed pion cm. angle 8 n = 90°. 
As shown by the solid curve in Fig. 4, our prediction for £(#„- = 90°) first tends to more 
negative values before bending over and becoming positive at large photon energies. It was 
found in Ref. |33| that in the threshold region this observable is very sensitive to the Mi_ 
multipole which strongly depends on the details of the low energy behavior of Roper res- 
onance, vector meson and FSI contributions. Therefore, a slight modification of one or all 
of these mechanisms can drastically change the photon asymmetry. As an illustration, our 
prediction for the energy dependence and angular distribution of S(6' 7r ) obtained with a 15% 
reduction of the Mi_ multipole, is shown by dashed curves in Fig. 4. This small modifica- 
tion of the low energy tail of the Roper resonance leads to positive photon asymmetries at 
all energies! 

In Table I, the ChPT predictions and the experimental values extracted from recent 
TAPS polarization measurements || are listed for comparison. Our predictions are in good 
overall agreement with the ChPT predictions |T(J and the TAPS results. However, there 
is a 15% — 20% difference in P 3 which leads to an underestimation of our result for the 
photon asymmetry, as shown in Fig. 4. Note that, in contrast to our model, P3 is essentially 
determined by a low energy constant in ChPT. 

Pion electroproduction provides us with information on the Q 2 = —k 2 dependence of the 
transverse E 0+ and longitudinal L 0+ multipoles in the threshold region. The "cusp" effects 

6 



in the Lq + multipole is taken into account in a similar way as in the case of Eq + , 
Re Lq+ = Re Z/g+ (IS) — a nN u c Re Lq+ (IS) 



\ 



£• (6) 



where all the multipoles are functions of total cm. energy E and virtual photon four- 
momentum squared Q 2 . It is known that at threshold, the Q 2 dependence is given mainly 
by the Born plus vector meson contributions in v!L, as described in Ref. [J30f . 



The validity of the approximation made in Eq. (||) is confirmed in Fig. 5 where the 
results of K-matrix approximation (solid lines) and full calculation (dash-dotted lines) agree 
with each other within a few percent. In Fig. 5 we also show our results for the cusp and 
FSI effects in the E 0+ and L 0+ multipoles for n° electroproduction at Q 2 = 0.1 (GeV/c) 2 , 
along with the results of the multipole analysis from NIKHEF || and Mainz || . Note that 
results of both groups were obtained using the P-wave predictions given by ChPT. However, 
there exist substantial differences between the P— wave predictions of ChPT and our model 
at finite Q 2 , as presented in Table II. To understand the consequence of these differences, 
we have made a new analysis of the Mainz data |J for the differential cross sections, using 
our DM prediction for the P-wave multipoles instead. The S"-wave multipoles extracted 
this way are also shown in Fig. 5 by solid circles and listed in Table III. We see that the 
results of such a new analysis gives a Eq + multipole closer to the NIKHEF data and in better 
agreement with our dynamical model prediction. However, the results of our new analysis 
for the longitudinal L 0+ multipole stay practically unchanged from the values found in the 
previous analyses. Note that the dynamical model prediction for L 0+ again agrees much 
better with the NIKHEF data. 

In Figs. 6 and 7, our model predictions (dashed curves) are compared with the Mainz 
experimental data || for the unpolarized cross sections da/dQ = dax/dQ + e da^/dVL, and 
for the longitudinal-transverse cross section daTi/dQ. Overall, the agreement is good. The 
solid curves are the results of our best fit at fixed energies (local fit) obtained by varying 
only the E Q+ and L 0+ multipoles. We have found that the differences between the solid and 
dashed curves in Figs. 6 and 7 are mostly due to the difference in the L 0+ multipole (see 



also Fig. 5). 

Finally we discuss a new version (hereafter called MAID) |53| of the unitary isobar model 
developed at Mainz (hereafter called MAID98) JJIJ, which is currently being intensively used 
for the analysis of pion photo- and electroproduction data. In this model the FSI effects 
are taken into account using the K-matrix approximation, namely without the inclusion of 
off-shell pion rescattering contributions (principal value integral) in Eq. (fj). As a result, 
the S-, P-, D- and F-waves of the background contributions are defined in MAID as 

t*(MAlD) = exp (z<y cos6 a v*(q E , k). (7) 

However, as we have found above, dynamical model calculations show that pion off-shell 
rescattering is very important at low pion energies. The prediction of MAID for tt° photo- 
production at threshold, represented by the dotted curves in Figs. 1 and 2, lies substantially 
below the data. It turns out that it is possible to improve MAID, in the case of ir° produc- 
tion at low energies, by introducing a phenomenological term and including the cusp effect 
of Eq. (|5|). In this extended version of MAID, we write the E 0+ (ti p) multipole as 

Re Ejf = Re E^f (MAID98) + E cusp (W, Q 2 ) + E corr (W, Q 2 ) , (8) 

where 



E cusp = -a nN u c ReE%Z (MAID98 



\ 



^ 2 



1-73- ( 9 ) 



UJ 2 



The phenomenological term E corr which emulates the pion off-shell rescattering corrections 
(or pion-loop contribution in ChPT) can be parameterized in the form 

E corr (W } Q 2 ) = {l + A B2ql)2 Fd(Q 2 ) , (10) 

where Fa is the standard nucleon dipole form factor. The parameters A and B are obtained 
by fitting to the low energy tt° photoproduction data: A = 2.01 x 10 _3 /m 7r + and B = 
0.71 fm. 

In summary, we have shown that within a meson-exchange dynamical model [EH 



one 



is able to describe photo- and electroproduction in the threshold region in good agreement 



with the data. The model has been demonstrated p4| , p5| to give a good description of most 
of the existing pion electromagnetic production data up to the second resonance region. The 
success of such a model at intermediate energies is perhaps not surprising since unitarity 
plays an important role there. However, it is not a priori clear that our model should 
also work well near threshold, even though we do start from an effective chiral Lagrangian. 
In principle, crossing symmetry is violated and the well-defined power counting scheme in 
ChPT is lost by rescatterings. As a matter of fact, previous similar attempts have failed 



22| . |24j| . It is easy to understand the difference between our present calculation and the 

are 



results of Ref. p2) by the fact that the off-shell behavior of the nN models used [13 



quite different from each other. The difference between our results and those of Ref. [24 



probably arises, in large part, from different off-shell prescriptions used for the transition 



potential t> 77r , because the meson-exchange nN model used in Ref. [24| is very similar to 



the one used in this study. On the other hand, meson-exchange models fll6| , |T7|l have also 
been shown to give a good description of low energy irN data, in addition to an excellent 
agreement with the data at higher energies. It is therefore assuring that similar success can 
also be achieved for the pion electromagnetic production. 

The largest discrepancy between our results and the data is in the P3 amplitude where 
our prediction underestimates the data by about 20%. As a consequence, our prediction for 
the photon asymmetry has the opposite sign as observed in the experiment ||. However, 
we have found that in the threshold region, the photon asymmetry is very sensitive to many 
ingredients of the theory, e.g., vector mesons, FSI and Roper resonance contributions, a fact 
that deserves further studies. 

Finally, we found that the effects of final state interaction in the threshold region and in 
the case of tt° production, are nearly saturated by the single rescattering term. Therefore, 
the existing one-loop calculations in ChPT could be a good approximation to threshold 
neutral pion production. 
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TABLES 





Born 


oj + p 


FSI 


res. 


tot. 


ChFT 


Exp. 


-^0+ 


-2.46 


0.17 


1.06 


0.07 


-1.16 


-1.16 


-1.33 ±0.11 


Pl 


9.12 


-0.35 


0.15 


0.38 


9.30 


9.14 


9.47 ±0.37 


p 2 


-8.91 


0.21 


-1.32 


-0.13 


-10.15 


-9.7 


-9.46 ±0.39 


P3 


0.18 


4.61 


3.36 


1.20 


9.35 


10.36 


11.48 ±0.41 



TABLE I. E 0+ (in units W~ 3 /m i: +) at threshold and Pi, P 2 and P 3 (in units 10~ 3 g/m 2 + ). 
Contributions of Born terms (Born), vector mesons ( tu + p), pion rescattering (FSI) and resonances 
(res.) are shown separately. The predictions of ChPT and recent experimental values are taken 
fromRef. iron and Ref. HI. 



Q 2 (GeV/c) 2 


0.0 


0.0 


0.05 


0.05 


0.1 


0.1 


Model 


DM 


ChPT 


DM 


ChPT 


DM 


ChPT 


-^0+ 


-0.92 


-0.96 


1.08 


0.45 


2.48 


1.60 


M 1+ 


1.11 


1.19 


1.70 


2.16 


1.93 


2.74 


Afi_ 


-0.59 


-0.49 


-0.96 


-0.72 


-1.07 


-0.74 


Pi+ 


-0.02 


-0.02 


-0.03 


-0.01 


-0.03 


-0.01 


Lq+ 


-2.71 


-1.61 


-1.74 


-1.52 


-1.13 


-1.31 


Li- 


-0.14 


-0.52 


-0.05 


-0.60 


-0.01 


-0.51 


Li+ 


-0.01 


-0.01 


-0.01 


-0.02 


-0.01 


-0.02 



TABLE II. Comparison of the S- and P-wave multipoles (in units 10 3 /?n, 7r +) for thep(7*,7r°)p 
reaction at Q 2 =0, 0.05 and 0.10 (GeV/c) 2 , obtained in our model and ChPT @ at AW= 2.5 

MeV. 
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AW 


Eq+ 


Lo+ 


X 2 /d.o.f. 


0.5 


2.28 ±0.36 


-1.34 ±0.06 


1.49 


Ref. § 


1.96 ±0.33 


-1.42 ±0.05 


1.29 


1.5 


2.43 ±0.21 


-1.38 ±0.04 


1.26 


Ref. § 


1.82 ±0.19 


-1.41 ±0.03 


1.19 


2.5 


2.98 ±0.20 


-1.38 ±0.04 


1.60 


Ref. | 


2.12 ±0.17 


-1.36 ±0.05 


1.68 


3.5 


2.61 ±0.22 


-1.39 ±0.03 


1.70 


Ref. § 


1.52 ±0.18 


-1.27 ±0.03 


1.84 



TABLE III. Values of the Eq + and Lq + (in units 10 3 /m 7r +) for the p(7*,vr°)p reaction at 
Q 2 =0.10 (GeV/c) 2 , obtained after the local fit to the differential cross sections measured in Ref. 
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FIGURES 
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FIG. 1. Real part of the Eq + multipole for 7p — > 7r°p. The dotted curve is the result 
obtained without FSI, the dash-dotted and short-dashed curves are the FSI contributions from 
charge exchange ir + n and elastic ir°p channel, respectively. The solid and long-dashed curves are 
the total results (see Eq. (0)) obtained with the full matrix t n ^ and with t n N replaced by v n N, 
respectively. 
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FIG. 2. Real (lower panel) and imaginary (upper panel) parts of the Eq + multipole for 
7P — > 7r°p. The dashed and solid curves are the full results obtained without and with isospin 
symmetry assumption, respectively. In the latter case, the pion mass difference effect is taken into 
account using Eq. (||). The dash-dotted curve is the result of ChPT [1C]. Data points from Ref. 
1(A), Ref. 0(.), and Ref. H(o). 
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FIG. 7. Same as in Fig. 6 for the transverse-longitudinal cross section daTi/dO,. 
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